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The three-compartmental demographic model by Korotaeyv-Malkov-Khaltourina, connecting 
population size, economic surplus, and educational level, is considered from the point of view of 
dynamical systems theory. It is shown that there exist two integrals of motion, which enable the 
system to be reduced to one non-linear ordinary differential equation. The study of its structure 
provides analytical criteria for the dominance ranges of the dynamics of Malthus and Kremer. 
Additionally, the particular ranges of parameters enable the derived general ordinary differential 
equations to be reduced to the models of Gompertz and Thoularis-Wallace. 
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I. INTRODUCTION 


The Population Division at the United Nations de¬ 
clared that “Timely and accurate information about pop¬ 
ulation trends continues to be in high demand” since 
it “is needed for the formulation and implementation of 

S ilicies and programmes in almost all areas of public life” 

. Correspondingly, the mathematical modeling of de¬ 
mographic processes dates back to the work of Malthus, 
Verhulst, Pearl, and Reed, and it remains a challenge for 
non-linear dynamics aimed at reproducing and explain- 
in g m odern population data and predicting future trends 

m- 

Two of the most notable features of global demogr^hy 
are the hyperbolic growth of the global population Q 


and the fast stabilizing balance between actual economic 
levels and population sizes. 

On the other hand, the demographic transition is a suf¬ 
ficiently more complex process. Its description requires 
taking into account several additional social factors lead¬ 
ing to a saturation of the population, and, correspond¬ 
ingly, more sophisticated models of non-linear dynamics. 
From the latter point of view. Ref. Q proved the exis¬ 
tence of the unified scaling approach not only to the gen¬ 
eralized model of von Foerster et al. [a(A^) = with 
arbitrary a] but also to the Thoularis-Wallace model Q 
of generalized logistic growth 


dN 

dt 


= bN° 



( 4 ) 


which is traced from I960-I970th back to prehistoric 
times [ 1 , 01 , and the global demographic transition, which 
slows down this population growth and tends to stabilize 
the overall number of people 0 , 0 ] . 

The hyperbolic function ([T|) is a solution of the equa¬ 
tion 


where K is the carrying capacity and a, /?, and 7 are the 
adjustable parameters. However, the developed mathe¬ 
matical framework leaves open the question of the origin 
of non-integer values of these parameters. 

Some recent sociological studies (lol - [ll| argue that ed¬ 
ucation level is one of the key factors that effectively 
stops world population growth. The primary effect is 
a lower overall fertility rate for well-educated women, 
which could dominate other economic and social factors. 

At the same time, Korotayev, Maikov and Khaltou- 
rina proposed [l3l| a relatively simple model that connects 
these three factors in a natural way: 


which was proposed in 0 considering the Malthusian 
dynamics 



( 3 ) 


with the population-dependent growth rate a = hN, 
which was hypothesized as an effect of people’s coop- 
erativity. A more precise background of this effect was 
revealed by Kremer 0, who argued that the functional 
form ([2]) originates from the Malthusian economic growth 
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aNS{l - L), 
bNS, 

cSL{l-L). 
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Here the population growth ([S]) satisfies Malthusian dy¬ 
namics with the growth’s rate supported by economic 
surplus S and restricted by education level L (the au¬ 
thors of m considered literacy to be such an indicator; 
however, some other educational determinants could be 
considered without loss of generality). The economic de¬ 
velopment (|ni) admits Kremer’s assumption, and educa¬ 
tion level satisfies the Verhulst form, in which the growth 
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rate is determined by the surplus. The coefficients a, b, 
and c are constants. 

It has been shown 0 that numerical solutions of the 
system ©-([T]) correspond to global statistical data with 
sufficiently high accuracy. However, all of these conclu¬ 
sions, including the detection of an interval, where N 
mimics hyperbolic growth (ED were the results of numer¬ 
ical simulations only. 

Thus, the main goal of the present work is to explore 
Eqs. (©-(IT]) analytically as a non-linear dynamical sys¬ 
tem. The principal questions are as follows: (i) What 
is its actual dimension, i.e., the existence of integrals 
of motion, that determines the degrees of freedom for 
an adjustment to the demographic data? (ii) Is it re¬ 
ally reducible to the classical models of demography and 
population dynamics in general, and what are the corre¬ 
sponding ranges of population size? 


II. RESULTS 


where 


So) = So- ^N 2 Fi (l, I + AiVo“) • (H) 

Therefore, both L and S can be substituted into 
Eq. (O as functions of population number, and the sys¬ 
tem Eq. (O by Eq. ED will be reduced to the unique 
equation with respect to this variable: 

+abfiN {1 - XN^) . (12) 

It should be supplied with an initial value N{0) = Nq. 

The hypergeometric function is given by the Gauss se- 
ries, 


2-^1 


^ Cl a , ,, 

1, 1 H—; XN 

^ c c 




r(i + f) “ r(f+ n) (Aivf)" 

r(f) h r(i + f + n) ’ 


A. Basic reduction of ordinary differential 
equations 


Since the system ©-([T]) contains similar products of 
the same variables, it is natural to examine scaling func¬ 
tions and integrals of motion considering derivatives of 
?V, S, and L with respect to each other, excluding time. 

It should be pointed out that one integral was found 
directly by the authors of the model dividing Eq. dSD by 
Eq. 0 that results in the power law dependence L = 


ATVa, where the factor A = LqNq “ is determined by the 
initial conditions. 

At the same time, this procedure can be successfully 
continued, i.e., the system ©-(ITl) can be reduced to one 
equation in an explicit form. Substitution of the revealed 
function L = L{N) into Eq. (I5|), 
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aNS{l-XN^), 

bNS, 


( 8 ) 

(9) 


which can be simplified using the functional relation 
r(l + z) = zT{z). 

As a result, Eq. (TT^ takes the form 


dN 

dt 


OO 

abfiN + bN^ + bN^ 

n—1 


{xn^Y 


(l-A7Vt), 


(13) 

which is quite demonstrable for a qualitative analysis. 

The second multiplier provides a delimiter, which de¬ 
termines the stable maximal population’s size Ns = 
^-i/(c/“) corresponding to the education level L = 1. 

The first multiplier consists of three terms. Each one 
corresponds to the basic demographic process: (i) abfiN 
is the simplest Malthusian growth; (ii) bN^ satisfies Kre- 
mer’s economic-demographic model; and (iii) the last se¬ 
ries describes the influence of education level on the rate 
of population growth. Each term dominates within its 
own range of N. 


B. Discussion of particular cases 


and dividing Eq. ([5|) by Eq. (O leads to the equation 

which allows for the separation of variables. 


N 


dN 


No 


1 - XNt 


= -^{S-So) 


The integral on the left-hand side has a closed-form ex¬ 
pression via the Gauss hypergeometric function 2 E 1 141. 
Thus, the solution explicitly representing A as a function 
of S, is 

S=-N2Fi (l,-;l + -;AiV^) +^i{No,So), (10) 

a \ c c / 


Small populations with an education level far from sat¬ 
uration {N « afi and N « grow exponen¬ 

tially as A oc exp{abfit) since Eq. (fOl) reduces to the 
simplest Malthus law, dN/dt = ab^N. 

The transition from exponential to hyperbolic growth 
occurs for N > a/x, but still N « X~^C{c/a)_ com¬ 
plex interplay of education level and population size is 
valuable at the area of the demographic transition only, 
when XN‘^/°' tends to unity. 

Note also that the global dynamics of the world’s pop¬ 
ulation exhibited purely hyperbolic behavior ii for the 
extremely long time interval before the end of the 20th 
century. This kind of solution is provided by Eq. E3D 
for XN^/°‘ if /X = 0. The last equality corresponds to the 
quite notable invariance property of the Cauchy problem 
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stated for the full system ([S])-©. This condition means 
that the equality 

S = -N2Fi( I,l) (14) 

a V c c / 

holds for all moments of time, from the initial one up to 
an arbitrary one. This follows from m and (EH) with 
the substituted value /r = 0 . 

The invariant equality (1141) corresponds explicitly to 
Kremer’s basic assumption that the population adjusts 
instantaneously to an economic indicator, but with the 
additional factor depending on education level L. This 
factor varies from 2 T 1 (l, f; 1 + f! O) =1 (and it is close 
to 1 for L « 1) up to 2F1 ( 1 , 1 + ^; 1 ) = a/c. 

Thus, the asymptotic balance between economic in¬ 
dicator S and population size N is represented by the 
values of the invariant El as 5 = 6a for the low- 
educated population (L = 0) and S = c~^N for the 
high-educated population [L = 1). Since the coefficient 
of proportionality is given as a ratio of economics and 
population growth rates h/a in the former case and as 
the inverse knowledge growth rate c~^ only in the lat¬ 
ter, the formula (ITTl) qualitatively expresses the typical 
characteristics of the transition from an industrial soci¬ 
ety to a post-industrial one [Il[i3. In addition, this 
effect, in conjunction with the valuable deviation of the 
Gauss hypergeometrical function from unity, should be 
taken into account in the solution of Eq. El), when JV 
tends to A“^. Thus, its behavior encompasses the coex¬ 
isting transitions in the demographic curve and in the 
socio-economical characteristics. 

Consider now some particular cases of parameters that 
reduce Eq. El ( or, equivalently, Eq. El) to known 
equations of population dynamics. 

The Gauss hypergeometric function has the following 
expression in elementary functions for a/c = 1 0: 

2 E 1 (1,1; 2; AN) = log (1 - AiV). 

Therefore, Eq. (IT3l) written with respect to dimensionless 
n = 1 — AN, takes the form 

flT) 1 — 71 

— = 6n [log[exp(a/r)] -b log(n)] . (15) 

Such an equation (the Gompertz equation multiplied by 
the Verhulst factor) was considered 0 in the context of 
the generalization of the Heumann-Hotzel model, which 
allows for a reproduction of the drosophila fly population. 

For n « 1, Eq. El reduces to the original Gompertz 
model 

= j [log[exp(a^)] -b log(n)] (16) 

with the solution 


where /c is a constant depending on a/i and initial condi¬ 
tions. 

Since it was proposed to discuss mortality rate 00, 
the description of the decay of the “illiterate” part of a 
population n is a natural association between models. In 
addition, the solution El demonstrates analytically a 
faster convergence to a stable stationary state in compar¬ 
ison with the logistic dynamics. This fact was mentioned 
as a result of numerical simulations for a/c = 1.32 in 
[l3l| . but it did not lead to any explanations or approx¬ 
imations there. Although the solution for a/c = 1.32 
differs quantitatively from the approximation a/c = 1 
El in a precise value, the principal qualitative behavior 
originates from the same hyperexponential dynamics. 

Another reduction to a known model can be evaluated 
for Eq. (fOl) taking a/c « 1 and /r = 0. One can neglect 
by the sum within the brackets there (in other words, to 
take 2 Ei(l,^;l-b|; AN o') Ri 1). Obviously, the remain¬ 
ing terms coincide with the Thoularis-Wallace model (|4]) 
with a = 2, P = c/a, and 7 = 1 . Thus, the question 
about the origin of an empiric fractional Verhulst fac¬ 
tor 0 can be reduced to consideration of growth rates 
in the extended system of ordinary differential equations 
for factors determining population growth. 


III. CONCLUSION 

The demographic dynamics, which satisfies Eqs. m- 
ED, is one-dimensional and can be reproduced by the 
single autonomous ordinary differential equation (El. Its 
coefficient can be taken as adjustable parameters, which 
can be determined from the known population trends. 
Toward that end, it is shown that the derived equation 
El reduces to the Malthus model or to the model of 
von Foerster et al. and Kremer when the population is 
far from the saturated state. 

The elimination of the Malthusian term is revealed 
as being directly connected with Kremer’s arguments on 
the balance between economic and population indicators 
with the addition of the people’s educational factor. The 
analysis of the latter allows for a demonstration of the 
transition from an industrial to a post-industrial society 
in terms of mathematical modeling. 

On the other hand, the interplay of these parameters 
with the growth rates of the original models and initial 
conditions places the described approach into the context 
of the methods of non-linear dynamics, which operate 
with the decomposition of complex differential equations 
of physico-chemical kinetics into a set of elementary poly¬ 
nomial ODEs [ 2 ^. 

In addition, the particular choice of parameters shows 
that the form Eq. EH is reducible to a variety of known 
non-trivial models of population dynamics, e.g., the mod¬ 
els of Gompertz and Thoularis-Wallace. Thus, it has 
potential applications beyond just the problems of de¬ 
mography. 
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